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The geometric relationships between ligands and the functional groups 
that bind ligands in soluble ligand-protein complexes have traditionally 
been deduced from distance constraints between pairs of NMR active 
nuclei spanning the ligand-protein interface. Frequently, the steep inverse 
distance dependence of the nuclear Overhauser effect (NOE), from which 
the distance constraints are derived, makes identification of sufficient 
numbers of constraints difficult. In these cases the ability to supplement 
NOE-derived information with distance-independent angular information 
can be very important Here, the observation of residual dipolar coup- 
lings from a-methyl mannose bound to mannose binding-protein in a 
dilute liquid crystalline medium has allowed the determination of a 
bound ligand's average orientation. The 3-fold rotational symmetry of 
mannose-b inding protein defines its orientational tensor and obviates the 
need to determine experimentally the protein's average orientation. 
Through superimposition of ligand and protein orientational tensors we 
describe the binding geometry of a-methyl mannose bound to mannose- 
binding protein. This new method is of general applicability to the study 
of ligands bound to proteins, and it is of particular interest when neither 
X-ray crystallography nor NOE techniques can provide sufficient infor- 
mation to describe binding geometries. 

© 1999 Academic Press 

Keywords: dipolar couplings; NMR; structure detennination; carbohydrate 
interactions; mannose-binding protein 



Introduction 

We illustrate here a new approach for the deter- 
mination of ligand geometry in protein binding 
sites that relies on orientational constraints derived 
from residual dipolar couplings observed in NMR 
spectra of partially ordered macromolecular com- 
plexes (Bothner-By, 1995; Prestegard, 1998; 
Prestegard et al., 1998; Tjandra & Bax, 1997; 
Tolman et al., 1995). The determination of the 
bound geometry of protein ligands is an important 
issue for improving our understanding of how 
structure relates to biological function, and also for 
building a basis for rational drug design. Many 
useful approaches to this determination have been 
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developed in the past, including X-ray diffraction 
of crystalline protein-ligand complexes and NMR- 
based transfer NOE studies of complexes in sol- 
ution (Ni, 1994). All the approaches are limited in 
certain ways. In the case of X-ray crystallography, 
complexes must be crystallizable; in the case of 
transferred nuclear Overhauser effect (NOE) stu- 
dies, ligands must be in rapid exchange. Transfer 
NOE studies, while giving sound information on 
ligand geometry, also suffer from a lack of infor- 
mation on the nature of ligand protein contacts. 
Information on these contacts can be provided by 
more conventional NOE studies using both protein 
and ligand resonances in some cases, but this infor- 
mation is often difficult to obtain because it 
requires assignment of protein as well as ligand 
resonances, and because it presumes sufficiently 
close approach of protein and ligand protons to 
allow efficient magnetization transfer (r<5 A). 

Cases where detennination of bound ligand geo- 
metry have been particularly difficult include those 
involving protein-oligosaccharide interactions. 
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3^$sETISiAte these processes is difficult oecause me 

ffe ^SgSSJS^ais are often part of the interface 
C^?ro2rir^d oUgosaccharide. Tne hydroxyl 
S^Ke^ theories in closest proximity to 
^^^oWbuT these exchange rapidly with 

water making their NMR reson- 
^hirf to observe. In tnese difficult cases, 
^ . , ^Jhdlpolar data can offer a valuable alterna- 
te^ ^wteelhe data can in pnnaple constrain both 
^CbSid^nd geometry and ligand onentahon 
5^3K? relative to the protein binding site. 
m-^^tJSaiLr target in the studies presented 
^CrSSding protein-A (MBP). MBP is 
a protein important in innate immune response. It 
binds sugars, including a-mannosides, that are fire- 
ouently exposed on the surface of invading patho- 
gens, and subsequently triggers a complement 
cascade (Dceda et al, 1987). Improper activation vm 
MBP has also been suggested as a factor in inflam- 
matory disease (Dwek, 1996). Crystal structures of 
several forms of MBP are available. The form we 
study here is a 149 residue truncated version that 
includes a caltium<ontaining carbohydrate recog- 
nition domain (CRD) and a short a-helical stem 
that initiates the formation of a homo-trimer by 
making a left-handed coiled-coil with helical stems 
from two other MBP monomers. The 3-fold sym- 
metry axis in the resulting homotrimeric structure 
is important for the interpretation of our results. 
The binding of several ligands to MBP and mutant 
MBPs that mimic selectin proteins has also been 
studied, both because of interest in MBP's ligand 
binding properties and those of proteins which 
show homology to the MBP CRD, such as the 
selectins. Crystallography has been the primary 
contributor in studies of MBP trimers, with struc- 
tures existing for bound galactoside and galactosa- 
mine, as well as sialyl Lewis-X and related ligands 
(Chang et al, 1994; Iobst et al, 1994; Kolatkar AR, 
1998; Ng & Weis, 1997; Weis & Drickamer, 1994). 
Additionally, there are structures of a shorter frag- 
ment of MBP which makes a homo-dimer; one of 
these includes oc-methyl mannoside in the carbo- 
hydrate binding site (Ng et al, 1996). The results of 
the above crystallography studies provide a useful 
basis for evaluating data from our new approach. 

The new approach that we present is based on 
the observation of residual dipolar coupling of 
directly bonded l H - l3 C pairs in a-methyl manno- 



side (AMM) when the mannoside is bound to MBP 
in a field-oriented aqueous liquid crystal. The use 
of residual dipolar coupling in the characterization 
of the structures of proteins and other macromol- 
ecules in solution has received a good deal of 
attention recently (Gore et al, 1998a; Tjandra & 
Bax, 1997) . When molecules are partially ordered 
in a magnetic field interactiorv vectors connecting 
magnetic nuclei in the molecules depart from a 
complete isotropic sampling of, afe cUrectkm fat 
space. This gives rise to a contribution to ipllttmg 
of NMR resonances that depend* 
properties of the interacting nuclei 
separation (r) and the angle (0) mat the vector 
makes with the magnetic field. When the nuclei 
are directly bonded as in a l H- u N ox l H- u C pain, 
the 1/r 3 distance-dependence of the interaction can 
be regarded as being determined by the bond 
length and the primary variable becomes the orien- 
tation of the bond relative to the magnetic field as 
characterized by the angle 9. The functional depen- 
dence of the dipolar couplings is (1/2(3 cos 2 
0 _ i)), where the parentheses denote an average 
over an orientation distribution, and the resulting 
interaction appears as an addition to the normal 
one bond scalar couplings. These through space 
interactions are always present; in fact, they are the 
basis of the common NOE interaction. However, 
the angular function above is averaged to zero 
when space is isotropically sampled and direct 
splittings of resonances only appear in partially 
oriented media. Fortunately, several means of pro- 
ducing field-ordered states are now available, 
including inherent orientation due to the large ani- 
sotropics in the magnetic susceptibilities of some 
molecules (Kung et al, 1995; Tolman et al, 1995), 
orientation due to the interaction of molecules with 
lipid bicelles that form field-oriented liquid crystals 
(Tjandra & Bax, 1997) and interaction with filamen- 
tous bacteriophage that also form field-oriented 
liquid crystals (Clore et al, 1998c; Hansen et al, 
1998). 

The interpretation of residual dipolar data has 
taken two routes, incorporation of individual bond 
constraints as penalty functions in simulated 
annealing protocols for molecular structure deter- 
mination (Clore et al, 1998a,b), and extraction of 
order tensors that directly describe the direction 
and level of ordering forces from the point of view 
of a coordinate frame fixed in a rigid molecular 
fragment (Losonczi et al, 1999; Losonczi & 
Prestegard, 1998b; Saupe, 1968). The latter 
approach which arose from substantial previous 
work in the liquid crystal field, proves to be par- 
ticularly valuable for the problem at hand. An 
order tensor is a 3 x 3 matrix with elements 
(1/2(3 cos 9 r cos 9j 8, 7 )> written in terms of direc- 
tion cosines that relate Cartesian axes of an arbitra- 
rily chosen molecular axis system to the magnetic 
field direction. Because the matrix is traceless and 
symmetric, there are only five independent 
elements. An order matrix in an arbitrarily chosen 
frame of a molecular fragment is difficult to inter- 
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pret, but if diagonalized, the five independent vari- 
ables turn into a principal order parameter, an 
asymmetry parameter and three Euler angles relat- 
ing the diagonal or principal frame, to the original 
molecular frame. If fragments of interest experience 
a common source of order, as they would if they 
were parts of a rigid ligand-protein complex, the 
directions and levels of the orienting force should 
appear the same from the point of view of each 
fragment when a model for the fragments is 
assembled with proper fragment orientations. This 
provides a test for the validity of proposed models 
for the orientation of a ligand (one fragment) in a 
protein binding site (another fragment). 

The key to using this test is having an approxi- 
mately rigid fragment with enough potential dipo- 
lar coupled pairs to determine the five elements of 
an order tensor. For a 15 N-labeled protein, pairs of 
'H- 15 N spins are clearly abundant and much of the 
protein backbone can be considered rigid. In oligo- 
saccharides, individual pyranose rings frequently 
prefer a *Ci configuration which can be considered 
to have approximately rigid geometry and ^H-^C 
pairs are abundant In o-methyl mannoside there 
are five directly bonded : H- l3 C pairs on the Hi! 
ring, two of which point in two unique directions 
and three of which depart slightly from a third 
direction. These pairs can provide the data for an 
order tensor determination. There are a number of 
ways of accurately detenruning the splittings of 
resonances from these pairs (Tolman & Prestegard, 
1996) but the easiest is simply to collect l H- I3 C het- 
eronuclear single quantum coherence (HSQC) spec- 



tra, without decoupling in one of the frequency 
domains. This is the approach we use here. We 
used uniformly ^XT-labeled ot-methyl mannoside in 
order to improve sensitivity. 

Results and Discussion 

Measurement of residual dipolar couplings in 
AMM complexed to MBP 

The l H- t3 C HSQC NMR spectra of l3 C-enriched 
AMM in the presence of an equimolar amount of 
MBP dissolved in a dilute bicelle medium are 
shown in Figures 1(a) and (b) at temperatures of 
25 °C and 39 °C, respectively. The data were 
acquired in the absence of ,3 C decoupling during 
acquisition, and the one bond ^-'H couplings are 
shown as frequency domain splittings in the pro- 
ton dimension. The magnitudes of these couplings 
and the corresponding errors extracted using a 
Bayesian parameter estimation method (An dree & 
Prestegard, 1998) are shown on the spectra. As one 
would expect for the case of isotropic tumbling at 
25 °C, all the couplings in Figure 1(a) axe similar 
and fall in the range of typical *Jch scalar coupling 
constants for non-anomeric sites in carbohydrates 
(140-150 Hz). On the other hand, couplings 
measured at 39 °C are more widely spread because 
of the presence of a residual dipolar contribution. 
Residual dipolar contributions for individual C-H 
bond vectors calculated from the differences in 
couplings measured at 25 °C (isotropic) and 39 °C 
(aligned) are shown in Table 1 (D^). 
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Table 1. Residual dipolar couplings for bound and free AMM 



Data Set 


C,-H, (Hz) 


C 2 -H 2 (Hz) 


C,-H, (Hz) 


C«-H« (Hz) 


C V H, (Hz) 


AMM + MBP + bicelle (D t *J 
AMM + bicelle (D„_) 
AMM-bound state (D*^) 


-10.8±0.7 
-3.5±0.7 
-22.6±2.2 


-15.2±0.7 
-7.3±0.7 
-28.7±2.2 


13.3±1.2 
9.2±0.9 
20±3.5 


16.6±l 
12.1 ±0.7 
23.9±2.9 


16.411.1 
7.8*0.5 
3O.5±3.0 



While residual dipolar U ^^ / J^ t, ^ t ^ 
are indeed observed in the AMM/MBP complex 
SLolved in a bicelle medium *ey need not ong- 
Sate fully from binding of AMM to a strongly 
Xned MBP molecule. We see only a single set o 
CTShces, but for a weak binding hgand in fast 
SdSr^that spends a portion of its time bound 
to^protein and the remainder free in solution, 
this will always be the case. Thus, observed split- 
tings will be a population weighted average of 
those induced in the free and bound states. 
Residual dipolar couplings have previously been 
observed infree carbohydrate molecules Q dissolved 
in bicelle media (Bolon & Prestegard 1998; Kiddle 
& Homans, 1998; Rundlof et al, 1998). While small, 
they can be significant. In order to quantify the 
contribution from the AMM-free state, the same 
experiments were repeated under identical con- 
ditions in the absence of MBP. The differences in 
measured couplings between 25 °C and 39 °C are 
also shown in Table 1 (D^). Again, couplings and 
associated errors were extracted using a Bayesian 
parameter estimation method. It is noteworthy that 
the magnitudes of couplings for various C-H bond 
vectors are generally smaller for free AMM 
samples, but the reduction varies from site to site. 
The different pattern of splittings assures us that 
we are not looking at a simple scaling of bicelle- 
induced order and that MBP has a significant 
direct effect on both extent and direction of AMM 
alignment. This direct binding effect was con- 
firmed in another set of experiments. Here, a molar 
ratio of 1:8 MBP: AMM was used, which would 
reduce the bound AMM state to approximately 
8%. Indeed, in this case where indirect effects 
would dorninate, the presence of MBP had little 
measurable effect on the observed residual dipolar 
couplings. 

Interpretation of residual dipolar contributions 
of a ligand in a protein-bound state will, in gener- 
al, require separation of the contribution from the 
i; free state using known binding properties. The 
r observed residual dipolar couplings measured in 
the AMM/MBP complex (Table 1, D^) will be a 
population weighted average from the free (N^) 
find bound (N^^) states, such that: 




^obs — NfreeDfree + NboundDbound 



(1) 



fraction of AMM free and bound (N^ and 
) can be calculated from the dissociation con- 
and known protein-ligand concentrations, 
a previously determined binding constant 
1 1 mM (Sayers, 1998), approximately 40% of 
is in fact bound to MBP leaving approxir 



mately 60% of AMM free in solution. Residual 
dipolar couplings originating from the bound state 
(Abound) calculated from equation (1) are shown In 
Table 1. Since these dipolar contributions originate 
from MBP's partial alignrne^ ^n>^ contain the 
desired orientational mfbrmation on AMM relative 
to MBP. This information is best extracted using an 
order matrix analysis. ' 



Order tensor calculation In AMM 

Using AMBER-rninimized structure coordinates 
for AMM (Woods et al, 1995) and five residual 
dipolar couplings measured in the sugar ring in 
simple bicelle solution, the five elements of the 
order matrix were determined using a singular 
value decomposition approach (Losonczi et al., 
1999). The initial molecular coordinate frame for 
AMM was defined by placing the y-axis along the 
CI -HI vector and the z-axis normal to the plane 
defined by Hl-Cl-OMe. Using this frame, the pos- 
ition of the principal order frame axes determined 
from residual dipolar couplings from free AMM 
(Dfrec Table 1) is illustrated in Figure 2(a) using a 
Sauson-Flaumsteed projection (Bugayevskiy & 
Snyder, 1995). The alignment tensor is extremely 
asymmetric and the directions of all three principal 
axes of the alignment tensor are well defined, 
except that inversion of any axis is allowed. The 
asymmetry is quantitatively defined by an asym- 
metry parameter eta (n) that can be written in 
terms of order parameters S ii (n = (S^ - S^J/SJ. 
It is equal to 0.8(±0.13). Asymmetric behavior is 
not unusual and has, in fact, been observed in 
many systems dissolved in bicelle media. Due to 
the high level of asymmetry, the direction of high- 
est order (SJ can alternate between two orthog- 
onal orientations (Figure 2(a)) for which 
assumes values similar in magnitude but opposite 
sign (S^ = + 0.0005 or -0.0005). Note that one of 
the possibilities for the higher level of order is 
nearly along the *-axis of the initial molecular 
coordinate frame (the center of the Sauson-Flaum- 
steed plot). Hence, it is nearly in the plane of the 
Hl-Cl-OMe fragment and perpendicular to the 
C-H vector. 

The five elements of the order tensor for AMM 
determined from couplings for the MBP-bound 
state (D^und) were also determined using a singu- 
lar value decomposition method, and the direc- 
tions of the principal axes are depicted in 
Figure 2(b). While the position of the z-axis, the 
most strongly ordered axis, is well defined, the 
position of the S xx and Syy axes are ambiguous. 



Coupling Constraints on Protein-Ligand Geometry 



111 



(a) 





Figure 2. (a) Sauson-Flaumsteed projection of the 
directions of highest order for oriented AMM without 
MBP. (b) Sauson-Flaumsteed projection of the directions 
of highest order for oriented AMM in the presence of 
MBP, with a dear depiction of axial symmetry. The 
molecular coordinate frame for AMM was defined by 
placing the y-axis along the Cl-Hl vector and the z-axis 
normal to the plane defined by Hl-Cl-OMe. 



This is a characteristic of an axially symmetric 
order tensor. The asymmetry parameters confirms 
this axially symmetric nature of the alignment ten- 
sor with T) = 0.24(±0.16). In addition, the level of 
order (S^ = - 0.0012(±0.00012)) is approximately 
2.5 times larger than that for AMM in free solution, 
and the direction of highest order is shifted 40° 
toward the molecular z-axis (out of the plane of to 
the Hl-Cl-OMe fragment). 

The direction of alignment for the bound ligand 
is most useful when referenced to the alignment of 
MBP. We know that the direction of highest order 
(S^) depicted in Figure 2(b) should coincide with 
the direction of highest order for MBP itself. Rotat- 
ing the molecular frame of the sugar to achieve 
coincidence then allows determination of the rela- 
tive orientation of AMM bound to MBP. In general, 
determination of the alignment from a protein is 
possible using l H- 15 N data as we have recently 
demonstrated for other systems (Fischer et ai, 



1999). However, careful consideration of MBP's 
structure, and in particular, its symmetry, allows 
us t omit this independent interpretation and 
indeed explain the appearance of axial symmetry 
depicted in Figure 2(b). 

MBP's alignment tensor a solution 
In symmetry 

The form of MBP used in these studies b a 
homo-trimer with s >fold symmetry axis and 
thre^equrvsknt tugl&tr£^^ by tte 
same symkneby operation as idepicted m Figure 3» 
The computed residual dipolar coupling ortrinat- 
ing from the AMM in the boimd state (l^ound) wiU 
thus be an average over three orientations that are 
related by MBrs 3-fold symmetry axis. As has 
long been recognized and used in studies of small 
molecules (Saupe, 1968) this averaging has a very 
specific effect on order matrices determined from 
averaged couplings; regardless of the nature of 
MBP's interaction with the bicelle medium, AMM's 
experimentally determined direction of highest 
order must always point along the symmetry axis 
of MBP and the resulting order tensor must always 
be axially symmetric (Al-Hashimi et al., 1999). 

In the absence of any large conformational flexi- 
bility and mobility in the bound state, the order 
parameters determined from bound AMM (S^ = 
- 0.0012(±0.00012)) will reflect MBP's order par- 
ameter. The fact that the principal order parameter 
is negative indicates that MBP's symmetry axis is, 
on average, perpendicular to the magnetic field 
and parallel with the normals of the bicelles in our 




Figure 3. Depiction of AMM in fast exchange with the 
three equivalent CRDs of MBP. The 3-fold axial sym- 
metry results in a direction of highest order along the 
rotor axis. Yellow spheres correspond to Ca 2+ , black 
and red shperes to carbon and oxygen, respectively, of 
AMM, and MBP is represented by ribbon diagram. 
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medium (these order with their normals perpen- 
dicular to the magnetic held; Sanders et al, 1994). 
This is as one would expect if interactions with the 
surfaces of bicelles dictated order of MBP. Axial 
symmetry should also propagate to the bound 
ligand and this is exactly what we observe in the 
order tensor determination using residual dipolar 
contributions from MBP bound AMM as shown in 
Figure 2(b). The direction of highest order, S^, in 
Figure 2(b) should now point along the symmetry 
axis of MBP. As we discuss below, the orientation 
of AMM deviates 





No crystal structure of trtmeric MBP-A with 
AMM in the binding site actually exists. However, 
AMM in the site of a dimer of the homologous 
MPB-C, missing the helical stem exists, and AMM 
in the trimer can be positioned by overlaying a 
carbohydrate-binding domain from the dimer on 
the three domains in the trimer structure (PDB 
entries 1KMB and 1RDL, rmsd 0.63-0.78 A). The 
orientation of AMM, when positioned so that the 
experimentally determined symmetry axes coincide 
with that of MBP cannot, however, be made to 
overlap with the orientation of AMM modeled in 
this way. It is displaced by approximately 40°. 
There are several possible explanations for this 
observed deviation. First, our use of the trimer as a 
model may be in error. Instead, an equilibrium 
among monomer, dimer and trimer forms may 
exist in solution. This possibility is unlikely, since 
line widths and diffusional properties of the con- 
struct used in previous studies are consistent with 
a trimer (Sayers, 1998). Moreover, the experimen- 
tally determined axially symmetric tensor is con- 
sistent with the expected 3-fold averaging. Second, 
the actual geometry of the trimer in solution may 
be different from what is depicted in the crystal 
structure. This is not out of the question; there are 
in fact X-ray structures of mutated MBPs in tri- 
meric aggregates which show substantial differ- 
ences in the relative placement of the CRDs, even 
though mutations are concentrated in loops at the 
extremities. For example, the helices near the tri- 
mer junction in PDB entries 1KMB and 1AFA are 
displaced from one another by more than a 1 A 
when the stems are aligned. Similarly, the presence 
of internal motion in each protein domain or in 
. AMM in the bound state can lead to deviations of 
i average orientations. Third, it may be that alternate 
binding geometries are possible! A second binding 
involving a second Ca 2+ , in each CRD has, in 
been observed in crystals containing a high 
1 of saccharide (Ng et al, 1996). If this site were 
a proper model would have to include 
of the orientation of the sugar in this 
mat the second and third cases would 
_ajdaUy symmetric nature of the deter- 
~£tfr but ^ direction of highest 



order will vary. A fourth possibility is that the 
binding geometry in the primary site can vary. 

In Figure 4 we show a geometry of AMM in the 
binding site of MBP (lKMB.pdb) that is consistent 
with our determined axis of symmetry. Using this 
refined structure, i-methyl mannose was docked 
to the binding calcium f the CRD of MBP. Orien- 
tational constraints required placing AMM's exper- 
imentally derived direction of highest order along 
MBPs 3-fold axis of symmetry. Subsequently, 
rotations around this' axis and translation along x, 
y, and v*fttk mamtairung a 16 A distance from 

AMM^MBP ^^mScttd tn Figure 4. 
Then fir* single violation of allowed van der 
Waab contacts (<Z5 A). This Is between HO-3 and 
Glul93 of MBP. The trimeric form used in our 
modeling has no ligand of AMM in the binding 
site and Glul93 coordination of the calcium com- 
pensates for the absence of sugar hydroxyl. Hence, 
it is quite likely that Glul93 moves on sugar bind- 
ing. MBP actually has a rather open binding site, 
in which coordination of the three and four oxygen 
atoms of the sugar to Ca 2+ are the primary inter- 
actions and considerable variation in the actual 
geometry of coordination may be possible. 

An independent determination of MBP's order 
tensor would remove some of the uncertainties in 
the above discussion. First, this would allow deter- 
mination of the relative domain orientations of 
MBP in solution. This can be achieved by equating 
the observed direction of highest order as seen 
from any single domain to the symmetry axis of 
trimeric MBP (Al-Hashirni et al, 1999). This pro- 
cedure, in fact, introduces a novel and simple route 
into determination of relative domain orientations 
in certain types of symmetry related homo-mul ti- 
mers such as MBP, a problem that so far has been 




Figure 4. Proposed mode of binding of AMM in MBP. 
The dark blue corresponds to Ca 2+ , and the light blue 
represents residues in the binding pocket of MBP. 
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very difficult to address using standard NMR tech- 
niques. Second, the determined direction of Highest 
order and the corresponding order parameters will 
have absorbed any Internal motional averaging 
contributions within MBP permitting direct com* 
partem with an ordtr tensor determined for AMM 
k tht bound sM04ditiorial averaging the 



diacctylated a-methyl rtumnoside. Flash chromatography 
followed by dcacrtyUHon with NaOMe/McOH pro- 
vided AMM in an 20% overall yield 



NMR spectroscopy 
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Material* "and'liWBhod* - : . 

Preparation of MBP-A 

MBP-A is an expression product from the pIN-mom- 
pA-2 plasmid (Drickamer, 1989), and was expressed as 
Jtescribed (Weis et d., 1991). Briefly, competent JA221 
cells were transformed with DNA plasmid and colonies 
selected from LB-amp plates for use in a starter culture 
in an LB-amp medium. Starter culture was grown to sat- 
uration overnight at 37 °C and used to inoculate five 
liters of LB-amp medium, which grew to ^550 of U.8. 
MBP-A production was then induced with ETG and 
additional growth allowed for 2.5 hours. Cells were har- 
vested by centrifugation at 4000 g, the pellet was resus- 
pended in 10 mM Tris-Cl (pH 7.0) and then lysed by 
sonication. The insoluble pellet from 18,000 g centrifu- 
gation was solubulized in 6 M guanidinium chloride, 
0.1 M Tris-Cl (pH 7.0), and clarified by centrifuging at 
138 000 g for 30 minutes at 4°C Dialysis of the centri- 
fuge extract with 25 mM Tris-Cl, 1.25 M NaCl, 25 mM 
CaCl 2 over 48 hours and centrifuging at 138,000 g pro- 
vided the dialysate which was purified by affinity chro- 
matography on a mannose-Sepharose column. Elution 
with 25 mM Tris-Cl, 1.25 M NaCl and 2.5 mM Na 2 EDTA 
followed by reconstitution by dialysis against 10 mM 
NaCl, 1 mM Tris-Cl, 25 mM CaCl 2 afforded 10 mg of 
MBP-A from a five liter growth culture. 

Preparation of AMM 

AMM was synthesized from [ l3 C 6 ] glucose by minor 
modification of a described method (Hare et al, 1993). 
Briefly, uniformly ,3 C-labeled D-glucose was protected as 
the penta-acetate with acetic anhydride, and the crude 
material brominated with HBr/AcOH. Reduction with 
Zn/aqueous AcOH provided the desired tri-0-acetyl-r> 
glucal, which was used in a Ferrier reaction with metha- 
nol/Dowex H + to give (a-methyl-2,3-dideoxy-4,6-di-0- 
acetyl-r>glucopyranoside. Purification by silica gel 
chromatography; followed by treatment of the anhyd- 
rous material with Os0 4 in pyridine afforded the desired 



hire wfth ^niaxfanal alignment 

trs were acquired using normal 

M ^ T (*h3*CI tingle quantum coherence 

ggS?S^ ^couplings 
JSSnisntfarSS direct proton dimension. Quadrature 
dScSnbi Ac r, evolution period was accomplished 
usinff gradient coherence selection. For all experiments, 
1024 points per scan, 16 scans per increment, 160 '1 incre- 
ments were acquired with a direct sweep width of 
3000 Hz and an indirect sweep width of 9000 Hz. Using 
the same acquisition parameters, another idenbcal Hb(£ 
experiment was acquired at a temperature of 25 C 
where the bicelle medium is in an isotropic solution 
state. 

Calculations of dipolar couplings and order tensors 

Dipolar couplings were calculated as the difference 
between the oriented couplings (Vch + D CH ) and the iso- 
tropic couplings (Vch)- m each case ' coup 1 "^ were 
extracted using a Bayesian time-domain NMR parameter 
estimation program Xrambo, using the method described 
(Andrec & Prestegard, 1998). This program 'is available 
on the Internet at http://tesla.ccrc.uga.edu. Typically, a 
2D HSQC data set with the desired splittings in the 
direct dimension was transformed and phased to yield a 
1024 x 256 real matrix. Single FID slices across the width 
of the resonance under investigation were analyzed inde- 
pendently. For every slice, the resulting frequency 
domain data were reverse Fourier transformed to gener- 
ate a ID time domain data set as input for Xrambo. The 
following model was used for the data. Each component 
of the doublet resulting from C-H couplings was given 
an identical linewidth and intensity, but an independent 
phase to circumvent any problems resulting from the 
presence of phase twist anomalies or dispersive contri- 
butions. Values for shifts, line widths, phases and intensi- 
ties were estimated and entered as starting parameters 
which were subsequently refined by Xrambo's Metropo- 
lis Monte Carlo method. This procedure was repeated 
for several indirect slices of the same doublet and the 
resulting rmsd across various slices was used as the 
uncertainty in measurement. The error analysis from this 
procedure was used to estimate the final precision of all 
measured residual dipolar couplings. The measured 
residual dipolar couplings and associated uncertainties 
along with AMM input coordinates obtained from an 
MD simulation were then used as input to a singular 
value decomposition program for the determination of 
order tensor elements (Losonczi et al, 1999). This 
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program is also available on the internet at http:// 
tesla.ccrc.uga.edu. 

Molecular dynamics stmulatl n 

MD simulations were performed with AMBER 
4.1 (Pearlman ei al, IW5) employing the alKatom 
GLYCAM93 parameter set for oligosaccharides (Woods 
it al., 1995) on a 12 processor SGI Origin 2000 comp- 
uter. AMM was solvated in an approximately 
32 A ic 32 A x 32 A box of 466 TIP3P water molecules 
and was minimized by steepest descent energy minimix- 
v v i#tofiarfn« a.diekctrk constant of IJO and a cut-off 

T?u&$7^ TlL± l— J~J «*ir SntoMwtinrMt r4 ft ft X Tnf+Ul 



utfajflrt were assigned at 5 K, and the simu- 



run, at 300 K at 1 atm. Selected torsions were 

> ;f extracted with me CARNAL module of AMBER 4.1 from 
— ^ trajectories to determine a closely related family of 
the moat populated structures. 

DocJdng ttudtoa 

First approximations of AMM's position in MBP were 
done wmS MIDAS Phis (Ferrin et al., 1988). Briefly, MBP 
with AMM modeled into the CRD were aligned with the 
helix parallel with the z-axis. A second AMM was 
iruitched and oriented according to the singular value 
decomposition analysis of the calculated dipolar coup- 
lings (a z - 10°, y + 40° rotation from the crystal struc- 
ture). This was then translated in x, y, and z directions 
while limiting rotations to the z-axis only. Solutions that 
maintained the three and four OH groups 2.6-2.8 A from 
the Ca +a responsible for binding and that rninimized van 
der Waals interaction with the protein were deemed 
plausible solutions. 
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